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Two microscopic models, UrQMD and QGSM, were employed to study the formation of locally 
equilibrated hot and dense nuclear matter in heavy-ion collisions at energies from 11. 6A GeV to 
160^1 GeV. Analysis was performed for the fixed central cubic cell of volume V = 125 fm 3 and for 
the expanding cell which followed the growth of the central area with uniformly distributed energy. 
To decide whether or not the equilibrium was reached, results of the microscopic calculations were 
compared to that of the statistical thermal model. Both dynamical models indicate that the state 
of kinetic, thermal and chemical equilibrium is nearly approached at any bombarding energy after a 
certain relaxation period. The higher the energy, the shorter the relaxation time. Equation of state 
has a simple linear dependence P — a(y / s)e, where a = c 2 s is the sound velocity squared. It varies 
from 0.12 ±0.01 at E lab = 11. 6 A GeV to 0.145 ±0.005 at E lab = 160,4 GeV. Change of the slope in 
a(y/s) behavior occurs at Ei ab = 40,4 GeV and can be assigned to the transition from baryon-rich 
to meson-dominated matter. The phase diagrams in the T — [ib plane show the presence of kinks 
along the lines of constant entropy per baryon. These kinks are linked to the inelastic (i.e. chemical) 
freeze-out in the system. 

PACS numbers: 25.75.-q, 24.10.Lx, 24.10.Pa, 64.30.-t 



I. INTRODUCTION 

Experiments on heavy-ion collisions carried out for 
the last two decades at GSI's Schwerionen Synchrotron 
(SIS), LBL's Bevalac, CERN's Super Proton Synchrotron 
(SPS), BNL's Alternating Gradient Synchrotron (AGS) 
and Relativistic Heavy Ion Collider (RHIC) have signifi- 
cantly helped us in understanding of properties of hot and 
dense nuclear matter. The collisions at top RHIC energy 
y/s = 200,4 GeV or at energy of coming soon CERN's 
Large Hadron Collider (LHC) ^/s = 5.5,4 TeV probe the 
domain of high temperatures and low net baryon den- 
sities, while the systems with lower temperatures but 
with much higher baryon densities should be produced 
in heavy-ion collisions at relatively moderate energies 
around E\ ab — 30^4 GeV accessible for future GSI's Facil- 
ity for Antiproton and Ion Research (FAIR) accelerator 
[l(. Most likely, the matter under such extreme condi- 
tions is composed of partons, i.e. quarks and gluons, in 
the phase of quark-gluon plasma (QGP), colored tubes of 
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chromo-electric field (or strings) , hadrons and their reso- 
nances. The question about the equation of state (EOS) 
of such substance remains still open. Present status of 
the nuclear phase diagram in terms of temperature T and 
baryon chemical potential is sketched in Fig. [TJ The 
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FIG. 1: (Color online) Present status of nuclear phase dia- 
gram in the T — /xb plane. 

highly anticipated transition between the hot hadron gas 
(HG) and the QGP is of first order for relatively dense 
baryonic substances only. With rising temperature and 
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dropping baryon density and chemical potential the tran- 
sition becomes of second order at the so-called tricritical 
point (TCP). After that it is transformed to a smooth 
crossover Q. Although the theory cannot localize the 
position of the TCP on T-fiB plane, lattice quantum chro- 
modynamic (LQCD) calculations indicate that it might 
be somewhere between the points with T w 160 MeV and 
Mb = 360 MeV [| or /i B = 470 MeV These values 
are close to the chemical freeze-out parameters obtained 
from the analysis of heavy-ion collisions at energies be- 
tween E lab = 11. 6 A GeV and E lab = AO A GeV within 
the statistical models [1, S 0| • They are close also to the 
temperatures and baryon chemical potentials in the cen- 
tral zone of heavy-ion reactions generated by microscopic 
transport models [H, [l(| EL E3 ■ Another interesting 
feature of the collisions at bombarding energies around 
40,4 GeV is the transition from baryon-dominated mat- 
ter to meson-dominated one. According to microscopic 
models, in gold-gold collisions at top AGS energy nearly 
70% of total energy is deposited in baryonic sector. At 
top SPS energy mesons are carrying 70% of the total 
energy, and at Ei ab w 40,4 GeV the energy parts of 
mesons and baryons are roughly the same. - The par- 
ticle composition is changing. Is it possible to trace con- 
sequences of this change in microscopic model analysis? 
To answer this question two transport Monte Carlo mod- 
els were employed: ultrarelativistic quantum molecular 
dynamics (UrQMD) model [l3| and quark-gluon string 
model (QGSM) The models use different mecha- 

nisms of string excitation and fragmentation. UrQMD 
relies on the longitudinal excitation, while the color ex- 
change scheme is employed in QGSM. Central gold-gold 
collisions with zero impact parameter b = fm were sim- 
ulated at bombarding energies Ei ab = 11.6,20,30,40,80 
and 160,4 GeV, respectively. Microscopic parameters 
related to quantities conserved in strong interactions, 
namely, the total energy, the net baryon charge, and 
the net strangeness extracted for a certain volume of the 
reaction were inserted into a system of nonlinear equa- 
tions to obtain temperature, baryon chemical potential 
and strangeness chemical potential of an ideal hadron gas 
in equilibrium. If the yields and transverse momentum 
spectra of particles obtained in microscopic simulations 
are close to that provided by the statistical model, the 
matter in the cell is considered to be in the vicinity of 
equilibrium. Then its equation of state and other ther- 
modynamic characteristics can be derived and studied. 

Relaxation of hot matter to equilibrium in the central 
cell of central heavy-ion collisions has been studied within 
the UrQMD model in @, i, [HI HO, d for energies rang- 
ing from 11.6,4 GeV at AGS to y/s = 200,4 GeV at RHIC, 
and, partially, within the QGSM [IS El- The size of the 
cell once chosen has been fixed throughout the system 
evolution. In the present paper we modify the analysis 
of the early stage in order to trace the expansion of an 
initially small area of homogeneity just after its forma- 
tion. The central volume was further sub-separated into 
smaller cells embedded within each other (" matryoshka- 



doll" structure). The transition of analysis from the 
smaller cell to the larger one was allowed if, and only 
if, the energy densities in both cells were the same. Re- 
gardless of the microscopic model applied for the actual 
calculations, the formation of (quasi)equilibrated state in 
the central cell at all bombarding energies in question is 
observed. The matter in the cell expands isentropically 
with constant entropy-per-baryon ratio. The isentropic 
regime arises even before the chemical and thermal equili- 
bration takes place. Due to coarse-graining of the central 
volume characteristic kinks in the temperature vs. bary- 
ochemical potential phase diagrams are found for both 
model simulations. This feature has not been seen in 
the previous studies because of the averaging of energy 
and baryon densities, in fact non-isotropically distributed 
within the relatively large volume. 

The paper is organized as follows. Similarities and dif- 
ferences of the microscopic string models chosen for the 
analysis are discussed in Sect. HU In Sect. IIIII criteria 
of thermal and chemical equilibrium are formulated, and 
Sect. HVI describes the statistical model of an ideal hadron 
gas used for the comparison with both microscopic mod- 
els. Section fVl presents the model study of the relaxation 
process in the cells with fixed and variable volumes. Spe- 
cial interest is paid to collisions at Ei ab — 20 A GeV and 
40,4 GeV, and obtained results are compared to that at 
neighbor energies. Conclusions are drawn in Sect. fVTl 



II. FEATURES OF URQMD AND QGSM 
A. Similarities of the microscopic models 

Both UrQMD and QGSM are formulated as Monte- 
Carlo event generators allowing to perform a careful anal- 
ysis of the measurable quantities by introducing all nec- 
essary experimental cuts. The models are designed to de- 
scribe hadronic, hadron-nucleus, and nuclear collisions in 
a broad energy range. In the hadronic sector both mod- 
els treat the production of new particles via formation 
and fragmentation of specific colored objects, strings. 
Strings are uniformly stretched, with constant string ten- 
sion k « 1 GeV/fm, between the quarks, diquarks and 
their antistates. The excited string is fragmenting into 
pieces via the Schwinger-like mechanism of qq and qq — qq 
pair production, and the produced hadrons are uniformly 
distributed in the rapidity space. 

To describe hadron-nucleus (hA) and nucleus-nucleus 
(A + A) collisions the momenta and positions of nucleons 
in the nuclei are generated according to the Fermi mo- 
mentum distribution and the Wood-Saxon density dis- 
tribution, respectively. The black disk approximation is 
adopted as criterion of interaction. It means that two 
hadrons can interact both clastically and inelastically 
if the distance d between them is smaller than \J a /tt, 
where a is the total cross section. Tables of the ex- 
perimentally available information, such as hadron cross 
sections, resonance widths and decay modes, are imple- 
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merited in the models. If this information is lacking, 
the one-boson exchange model, detailed balance consid- 
erations and isospin symmetry conditions are employed. 
The propagation of particles is governed by Hamilton 
equation of motion, and both models use the concept of 
hadronic cascade for the description of hA and A + A 
interactions. Note that such a rescattering procedure is 
very important in the case of relativistic heavy-ion colli- 
sions and is necessary for the thermalization of the fire- 
ball. Due to the uncertainty principle newly produced 
particles can interact further only after a certain for- 
mation time. However, hadrons containing the valence 
quarks can interact immediately with the reduced cross 
section a = u q N- The Pauli principle is taken into ac- 
count via the blocking of the final state, if the outgoing 
phase space is occupied. The Bose enhancement effects 
are not implemented yet. 



B. Differences between the microscopic models 



that the basic underlying principles and designs of the 
models are quite far from each other. By using both the 
UrQMD and QGSM for studies of the relaxation process 
in a broad energy range one can expect that the model- 
dependent effects, caused by application of a particular 
event generator, will be significantly reduced. 



III. STATISTICAL MODEL OF AN IDEAL 
HADRON GAS 

For our analysis of the thermodynamic conditions in 
the cell we use a conventional statistical model (SM) of 
an ideal hadron gas formulated in pioneering works of 
Fermi [22| and Landau [23[ . The statistical approach was 
successfully applied to the description of particle produc- 
tion in heavy-ion collisions from AGS to RHIC energies 
(see @ and references therein) . In chemical and thermal 
equilibrium the distribution functions of hadron species 
i at temperature T read (in units ofc = fcs = fi = l) 



The differences between the models for hadronic inter- 
actions arise on three stages. First stage is the formation 
of strings. The UrQMD belongs to group of models based 
on classical FRITIOF model [17], while the QGSM uses 
the Gribov Reggeon field theory (RFT) HEEl. In thc 
FRITIOF model the longitudinal excitation of strings is 
employed, and the string masses arise from momentum 
transfer. In the Gribov-Regge models the string masses 
appear due to the color exchange mechanism, and strings 
are stretching between the constituents belonging to dif- 
ferent hadrons. Longitudinal excitation of strings is also 
possible in the QGSM. This mechanism describes the 
processes of single and double diffraction. The second 
stage concerns the string fragmentation. The Lund JET- 
SET routine [2(|, used in the UrQMD, assumes that the 
string always breaks into a sub-string and a particle on a 
mass shell. In the QGSM the Field-Feynman algorithm 
[2ll ] with independent jets is applied. Therefore, the frag- 
mentation functions which determine the energy, momen- 
tum, and the type of the hadrons produced during the 
string decay, are different in the models. Thc third stage 
deals with the number and type of the stings produced in 
the collision. Due to the different mechanisms of string 
excitation and fragmentation, these numbers are also dif- 
ferent for two microscopic models in question. Last but 
not least, both models do not use the same tables of 
hadrons, chosen as discrete degrees of freedom. Whereas 
the UrQMD contains 55 baryon and 32 meson states to- 
gether with their antistates, the QGSM takes into ac- 
count octet and decuplet baryons, and nonets of vector 
and pseudoscalar mesons, as well as their antiparticles. 
Further details can be found in d3j and [14| . Recently, 
the QGSM has been extended by the implementation of 
a parton recombination mechanism (l6j . Since parton 
recombination plays a minor role for nuclear collisions 
at intermediate energies, the whole analysis of the relax- 
ation process is done for the standard QGSM. We see 



f(p,m l ) = 



exp 



T 



± 1 



(1) 



where p, mi, ti — \/p 2 + mf , and /i, are the full momen- 
tum, mass, energy, and the total chemical potential of the 
hadron, respectively. The "+" sign is for fermions and 
the "— " sign for bosons. Since in equilibrium the chemi- 
cal potentials associated to nonconserved charges vanish, 
the total chemical potential assigned to the i-th hadron 
is a linear combination of its baryon chemical potential 
/zb and strangeness chemical potential 



(2) 



with Bi and Si being the baryon charge and the 
strangeness of the particle, respectively. The isospin 
chemical potential (or, alternatively, chemical potential 
associated with electric charge) is usually an order of 
magnitude weaker than ^ B and Therefore, the depen- 
dence on this potential is disregarded in Eq. ([2]). Then, 
particle number density rij and energy density £j are sim- 
ply moments of the distribution function 



9, 



(2tt) 3 
9i 



f(p, mi)d 3 p , 



(3) 



^3 / ^p 2 + m? i f(p,m i )<Pp, (4) 



(2tt) e 



with gi being the spin- isospin degeneracy factor of hadron 
i. The partial hadron pressure given by the statistical 
model reads 



P = 



(^/ a^+V ^^- (5) 



The integrals in Eqs. ([3]) — © can be calculated numer- 
ically. Another way is to use a series expansion of Eq. (fTJ 
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in the form 22 



/(p,m 2 ) = £(Tl)" +1 exp 



Ei — p t 



, (6) 



which is inserted into Eqs. ©—([SI)- After some straight- 
forward calculations one gets 
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where ifi and -ft^ are modified Hankel functions of 
first and second order, respectively. The first terms 
in Eqs. (fZj) — (J9j) correspond to the case of Maxwell- 
Boltzmann statistics, which neglects the ±1 term in par- 
ticle distribution function (HI) - 

The entropy density in the cell is represented by a sum 
over all particles of the product f(p, rut) [1 — In /(p, m^)] 
integrated over all possible momentum states 



= - E £2 L /(p. ™») t ln /(p. m *) - !] p 2d p- ( 10 ) 



According to the presented formalism, the hadron com- 
position and energy spectra in equilibrium are deter- 
mined by just three parameters, namely, the tempera- 
ture, the baryon chemical potential, and the strangeness 
chemical potential. In order to define values of T, /ie, 
and /is one has to obtain the total energy density e, 
baryon density pe and strangeness density ps for a given 
volume from microscopic model calculations, and insert 
them as input parameters into the system of nonlinear 
equations 



Pb = ^2Bini(T,HB,Hs) , 

i 

Ps = y^,Sim(T,[i B ,ns) , 

i 

£ = E e '( T '^ B '^ s ) ' 



(11) 

(12) 
(13) 



where n,(T, fi B , fi s ) and £{(T, pb,Ps) are given by 
Eqs. ©-(H]). Since the particle data tables implemented 
in the microscopic models contain different numbers of 
hadrons, two versions of the SM with properly adjusted 
lists of hadron species are used, i.e. the number of 
hadronic degrees of freedom in the macroscopic model 
should correspond to that in the microscopic model. To 
decide whether or not the equilibrium is reached the cri- 
teria of the equilibrated state for open systems, discussed 
in the next section, should be applied. 



IV. CRITERIA OF THERMAL AND 
CHEMICAL EQUILIBRIUM 

Criteria of local equilibrium for open systems were for- 
mulated in 0], and we recall them briefly. Compared to 
a nonequilibrium state, the equilibrium is characterized 
by the absence of collective effects, like flow of matter 
or flow of energy. The fireball produced in heavy-ion 
collisions is always expanding both radially and longitu- 
dinally. Therefore, the centrally placed symmetric cell is 
chosen to diminish effects caused by nonzero collective ve- 
locity of any asymmetric or asymmetrically located cell. 
The cell should be neither too small to allow for the sta- 
tistical treatment, nor too large, - otherwise the homoge- 
neous distribution of matter may not be reached. Previ- 
ous studies d, H, [l(| HH, EH found that the cubic cell of 
volume V — 125 fm 3 centered around the center-of-mass 
of colliding gold-gold or lead-lead nuclei is well suited for 
such an analysis. Clearly, the relaxation to local equilib- 
rium cannot occur earlier than at a certain time needed 
for the Lorentz contracted nuclei to pass through each 
other and leave the cell 



„„ 2i? Az 
t eq > — H 

" 7/3 2/3 



(14) 



Here R is the nuclear radius, Az is the cell length in lon- 
gitudinal direction, (3 is the velocity of nuclei in the cm. 
frame, and 7 = (1 — . Quite unexpectedly, the 

reduction of the longitudinal size of the cell from 5 fm to 
1 fm does not automatically imply a faster equilibration 
in the smaller cell: the transition times are practically 
the same [l(| • This means that the transition to equilib- 
rium takes place simultaneously within a relatively large 
volume along the beam axis. 

Isotropy of the pressure gradients is a necessary con- 
dition for kinetic equilibration. Diagonal elements of the 
pressure tensor P{ x ,y, z } are calculated from the virial the- 
orem [H 



1 

W 



2 



(15) 



where V, rrii and pi are the volume of the cell, the mass 
and the momentum of the ith hadron, respectively. Fig- 
ure depicts the convergence of the transverse pressure 
in the cell to the longitudinal one in the UrQMD and the 
QGSM calculations. Both models claim that the pressure 
becomes isotropic at t < 10 fm/c after beginning of the 
collision. The time of convergence decreases with rising 
bombarding energy. The pressure calculated according to 
the statistical model is plotted onto the results of micro- 
scopic simulations also. The agreement between micro- 
scopic and macroscopic calculations is good for a period 
of about t = 8 — 10 fm/c. Then the matter in the cell 
becomes quite dilute, and the collision rate is not suffi- 
ciently high to maintain equilibrium anymore. However, 
the isotropy of pressure can be obtained, for instance, in 
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FIG. 2: (Color online) The longitudinal (3P Z , dashed curves) 
and the transverse (3P X , dash-dotted curves) diagonal com- 
ponents of the microscopic pressure tensor in the central 125 
fm 3 cell in (a) UrQMD and (b) QGSM calculations of central 
Au+Au collisions at energies from 11.6j4 GeV to 158A GeV. 
Asterisks indicate the pressure given by the statistical model 
and solid lines show the total microscopic pressure. 



a spherically expanding system of non-interacting parti- 
cles. To exclude such a situation from the analysis one 
has to impose two additional criteria concerning thermal 
and chemical equilibrium. 

For a closed system in equilibrium the distribution 
functions of particles are given by Eq. ([T]) with an unique 
temperature, so the hadron composition and energy spec- 
tra are fixed. In open systems neither the energy density 
nor the number of particles is conserved. Therefore, the 
snapshots of hadron abundances and energy spectra ob- 
tained at a certain time t should be compared with those 
corresponding to an ideal gas in equilibrium. The techni- 
cal procedure is simple. At the very beginning, the pres- 
sure gradients in transverse and longitudinal directions 
are considered. If the pressure isotropy is restored, say, 
within 10%-limit of accuracy, the densities of conserved 
quantities, i.e. energy, baryon charge, and strangeness, 
determined microscopically, (i) should be used as an in- 
put to Eqs. (|ll[) - (fT3"|) . The solution of this system of 
equations (ii) provides us with values of the tempera- 
ture, baryon chemical potential, and strangeness chemi- 
cal potential which fully determine the composition and 
spectra of particles. By (hi) a comparison of microscopic 



and macroscopic yields of the most abundant hadronic 
species one can decide whether or not the chemical equi- 
librium occurs, whereas (iv) the energy spectra of these 
hadrons should possess a common slope corresponding to 
1/T (thermal equilibrium). The similarity of the particle 
distributions means that our system is in the vicinity of 
equilibrium. At each subsequent time step the procedure 
described by (i)-(iv) is repeated. 



V. RELAXATION TO EQUILIBRIUM. RESULTS 
AND DISCUSSION 

A. Yields and energy spectra 

The yields of main hadron species, i.e. N, A, A + 
E, 7T, K and K in the central cell are shown in Fig. [3] 
for central Au+Au collisions at Ei a b = AO A GeV. For all 
particles, except pions, the agreement between the mi- 
croscopic and macroscopic estimates at t > 9 fm/c is 
good. Compared to the microscopic models, the number 
of pions is underestimated in the SM. The pion excess 
comes from the many-body decaying resonances, such 
as N*,A*,A*,lu, etc, and strings. After t = 10 - 13 
fm/c the many-body processes play just a minor role, and 
the pion multiplicity slowly converges to the equilibrium 
value. It looks like all species of the hadronic cocktail, 
except pions, are not far from the chemical equilibrium. 
It is well-known that the pure statistical model of an 
ideal hadron gas, which does not include effective chem- 
ical potential for pions or weak decays, systematically 
underestimates the pion yields compared to experimen- 
tal data. Nevertheless, the excess of pions in the model 
with short table of resonances, QGSM, is quite signifi- 
cant. This circumstance should affect the thermal spec- 
tra of all hadrons, provided the thcrmalization is reached. 

To verify how good the temperature is reproduced, the 
energy spectra dN/AnpEdE are displayed in Fig. 2J The 
Boltzmann fit to particle distributions is performed, and 
the SM calculations are plotted onto the microscopic re- 
sults also. Both in UrQMD and in QGSM the energy 
spectra agree well with the exponential form of the Boltz- 
mann distributions. Despite the good quality of the fit, 
the abundance of pions in particle spectrum leads to sig- 
nificant reduction of the effective temperature of the sys- 
tem within the QGSM calculations. Analytical estimates 
of the temperature drop (see Appendix\E$ are close to the 
temperatures extracted from the fit. It would be possible 
to diminish the pion yield by taking into account larger 
part of the resonance states, but our intention is to check 
the principal occurrence of the (quasi) equilibrium states 
in different microscopic models and to define the limits 
imposed on the effective equation of state. Note also, 
that significant part of the pion spectrum seems to be 
softer compared to other hadronic species. These pions 
are coming mainly from the decays of resonances and 
experience too few elastic collisions, that are necessary 
for their thermalization. Since the hadronic matter in 
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FIG. 3: Evolution of yields of hadron species in the central 
cell of volume V = 125 fm 3 in (a) UrQMD and (b) QGSM 
calculations (histograms) of central Au+Au collisions at 40^4 
GeV. Asterisks denote the results of the statistical model. 



the central cell reaches the state of thermal equilibrium, 
one can apply the mathematical apparatus formulated in 
Sect. IIHI and, finally, obtain the anticipated EOS. 



B. Evolution of the cell characteristics 

According to the information provided by Figs.[5]-|3J 
the appropriate time to start the study of thermodynamic 
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FIG. 4: (Color online) Energy spectra of N (Y), A (□), n 
(•), K (A), K (■), and A (O) in the central 125 fm 3 cell in 
(a) UrQMD and (b) QGSM calculations of central Au+Au 
collisions at 40^4 GeV at t = 13fm/c and t = lOfm/c, re- 
spectively. Lines show the results of the fit to Boltzmann 
distribution. 



conditions in the cell is t = llfm/c for the reactions at 
E lab = 20 A GeV and t = 9fm/c for E lab = 40A GeV. 
The input parameters obtained in the microscopic model 
analysis are listed in Tables HTHll together with the out- 
put thermodynamic characteristics given by the SM. Be- 
cause of the different number of hadronic states employed 
by QGSM and UrQMD, the tables of available hadronic 
degrees of freedom in the statistical model are adjusted 
properly. The only objects not taken into account in the 
SM are strings. The detailed analysis done in [25|, Hi| 
shows that string processes play only a minor role at such 
late times in the central part of the reaction. Less than 
5% of the total amount of hadronic collisions result to 
formation of strings. The strings produced at late time 
stages are quite light, and usually just one extra-particle, 
most commonly a pion, is produced after the string frag- 
mentation. This circumstance, however, may account for 
the pion overproduction (see Fig. [3]), since the inverse re- 
actions such as 3(or more)— >2 are not incorporated in the 
employed versions of both microscopic models. 

For both energies the baryon density in the cell at the 
beginning of the equilibrium phase is about 30% larger 
than the normal baryon density pg = 0.16fm~ 3 in the 
UrQMD calculations. Whereas QGSM allows for the 
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TABLE I: The time evolution of the thermodynamic char- 
acteristics of hadronic matter in the central cell of volume 
V = 125 fm 3 in central Au+Au collisions at bombarding en- 
ergy 20^4 GeV. The temperature, T, baryochemical potential, 
/Ub, strange chemical potential, ps, pressure, P, entropy den- 
sity, s, and entropy density per baryon, s/pB, are extracted 
from the statistical model of ideal hadron gas, using the mi- 
croscopically evaluated energy density, e cc11 , baryonic density, 
p B el \ an d strangeness density, p| e11 , as input. Of each pair of 
numbers, the upper one corresponds to the UrQMD calcula- 
tions, and the lower one to the QGSM calculations. 
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TABLE II: The same as Table U but for 40,4 GeV. 
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90.5 
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69.9 
73.8 
55.0 
59.0 
43.3 
47.8 



0.226 
0.290 
0.175 
0.219 
0.135 
0.170 
0.104 
0.130 
0.081 
0.101 
0.064 
0.077 
0.050 
0.060 
0.040 
0.047 
0.032 
0.038 
0.026 
0.031 
0.021 
0.026 



-0.0181 
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-0.0145 
-0.0041 
■0.0113 
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160.2 
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146.8 
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128.0 
126.3 
122.3 
120.5 
117.2 
115.2 
112.0 
110.2 
107.0 
105.7 
102.4 
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341.6 
401.9 
354.2 
417.9 
363.7 
434.5 
374.7 
447.2 
390.1 
460.3 
404.0 
465.9 
419.0 
471.9 
426.6 
479.7 
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489.8 
457.3 
499.7 
469.8 
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75.5 
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71.6 
100.4 
67.4 
94.7 
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79.4 
50.8 
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46.3 
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67.2 
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34.2 
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4.23 
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3.25 
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2.53 
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1.90 
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1.15 
0.93 
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0.74 
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18.69 
15.01 
18.60 
14.85 
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14.61 
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14.53 
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14.85 
18.43 
15.16 
18.81 
15.38 
18.69 
15.39 
18.40 
15.48 
18.34 
15.31 



production of hot equilibrated matter with a density of 
Pb = 1-8/90, much higher nuclear densities obtained in 
microscopic simulations have been reported [27] ■ One 
has to bear in mind two important things concerning such 
density estimates. Firstly, they are very sensitive to the 
volume of the test-system, especially at the initial stage 
of the collision. As seen in Fig. H^a) and Fig. Ufa), the 
baryon density in both models cannot exceed 5 po in the 
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FIG. 5: (Color online) The total energy density e versus 
baryon density pb obtained in the central UrQMD cell of 
volume (a) V = 125 fm 3 and (b) V = 0.125 fm 3 during the 
time evolution of central Au+Au collisions at energies from 
11.6.A GeV to 158A GeV. Dashed lines correspond to the 
non-equilibrium stage of the reaction, solid lines represent the 
equilibrium phase. 



central cubic cell with volume V — 5 x 5 x 5 fm 3 regardless 
of the bombarding energy, while for the smaller cell with 
volume V sma ii = 0.5x0.5x0.5 fm 3 the baryon density can 
be as high as 20 po in the calculations within the same mi- 
croscopic models, see Fig. [Sfb) and Fig. [U^b). Secondly, 
such high values of the ps should be treated with great 
care. The accelerated cold nuclear matter is automati- 
cally "compressed" in the calculations by the 7- factor. At 
the initial stage of a nuclear collision one deals with two 
opposite fluxes of Lorentz-contracted nucleons which just 
start to interact with their counterparts. Although the 
calculated baryon densities are huge, especially for the 
small cell, this is a purely kinematic effect, since the sys- 
tem is far from local equilibrium. The numbers become 
meaningful only when the equilibration takes place. 

Another interesting effect is the negative (though 
small) net strangeness density in the cell throughout the 
evolution of the system depicted in Fig. [7] The result is 
pretty insensitive to the size of test-volume and can be 
explained as follows. Strange particles are always pro- 
duced in pairs, so the net ps is zero. At energies about 
40,4 GeV and below kaons emerge predominantly with 
lambdas and antikaons. Because of the s-quark in its 
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FIG. 6: (Color online) The same as Fig. [5] but for QGSM 
calculations. 



FIG. 7: (Color online) Time evolution of net strangeness 
density in the central 125 fm 3 cell in (a) UrQMD and (b) 
QGSM calculations of central Au+Au collisions at energies 
from 11. 6A GeV to 158A GeV. Lines are drawn to guide the 
eye. 



composition, kaons have significantly smaller interaction 
cross-section with baryons at p < 2GeV/c compared to 
antikaons, which carry the s-quark. Therefore, K leave 
the central cell with positive net baryon charge easier 
than A or K thus resulting to negative net strangeness. 
At RHIC energies the B — B asymmetry in the cell is 
much less pronounced, and the net ps is very close to 
zero @, EH- 

Here we distinctly see the role of hadronic degrees of 
freedom. Despite the net baryon density is about 15% 
larger in the QGSM calculations than in the UrQMD 
ones, the absolute value of the net strangeness density is 
almost 30% higher in the UrQMD cell as compared to 
that in the QGSM. Extra-strangeness is deposited in the 
resonance sector, mainly in A* and K* . Although the net 
ps i n the cell shown in Fig. [7] quickly drops almost to zero 
after t = 6fm/c, its relaxation proceeds slower than that 
of the net baryon density. Figure [8] displays the instant 
rise of the ratio f s = —ps/pB with time t attributed 
to both microscopic models. The role of the small non- 
zero net strangeness is not negligible. The difference in 
particle spectra and, especially, in particle ratios can be 
about 15% if one performs the SM calculations with 
essentially zero net strangeness. 



C. EOS in the cell 

Isentropic expansion of relativistic fluid is one of the 
main postulates of Landau hydrodynamic theory [23| of 
multiparticle production. We cannot prove or disprove 
this assumption in microscopic simulations for the whole 
system, simply because a global equilibrium is not at- 
tained. Though conditions in the cell are instantly chang- 
ing, it is possible to check the behavior of the entropy per 
baryon. Within the 5% accuracy limit, this ratio is nearly 
conserved in the equilibrium phase of the expansion, see 
Fig. [5J The entropy densities obtained for the cell in both 
models are very close to each other, but, because of the 
difference in net-baryon sector, the ratio s/ pe in UrQMD 
is about 15-20% larger than that in QGSM. Together 
with the pressure isotropy, the conservation of entropy 
per baryon supports the application of hydrodynamics. 

Any hydrodynamic model relies on the equation of 
state, which links the pressure of the system to its energy 
density. Otherwise, the system of hydrodynamic equa- 
tions is incomplete. The corresponding plot with micro- 
scopic pressures P m ic{£) is presented in Fig. [Till whereas 
the macroscopic pressures obtained from the SM fit are 
shown in Fig. [TJJ In the last plot the dependence of 
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FIG. 8: (Color online) The same as Fig.[7]but for strangeness 
per baryon, f s = —ps/ps- Lines are drawn to guide the eye. 

pressure on energy density is remarkably linear for both 
models for all energies in question. Thus the EOS has a 
rather simple form 

P(e) = c 2 s e , (16) 

where the sonic velocity in the medium c s — (dP/de) 1 / 2 
is fully determined by the slopes of the distributions 
P(e). However, if the pressure is determined microscop- 
ically and not via the distribution function, the falloff of 
pressure with decreasing energy density proceeds slightly 
nonlinearly. This feature can be seen distinctly for top 
SPS energy in the QGSM calculations. Therefore, for 
both models we averaged the slopes of the P vs. e distri- 
butions over the whole period of the equilibrated phase 
(see Fig. [10]). It should be noted that due to the aver- 
aging over time, respectively energy density, the values 
do not represent the maximal values for c 2 which are 
reached in the corresponding reactions. They are actu- 
ally lower, since also energy densities below the critical 
energy density of about 0.8GeV/fm 3 contribute to the 
average. 

The extracted values of the c 2 are presented in Fig. [T^J 
For the UrQMD calculations the velocity of sound in- 
creases from 0.13 at E lab = 11. 6 A GeV to 0.146 at 
Ei a b = 158, A GeV, and saturates at c 2 = 0.15 for RHIC 
energies, y/s = 130 A GeV and y/s = 200 A GeV [HI- In 
QGSM calculations the averaged sound velocity is about 



FIG. 9: (Color online) Time evolution of entropy per baryon 
S/ps in the central 125 fm 3 cell in (a) UrQMD and (b) 
QGSM calculations of central Au+Au collisions at energies 
from 11.6j4 GeV to 158 A GeV. Dashed lines correspond to 
the non-equilibrium stage of the reaction, solid lines repre- 
sent the equilibrium phase. 

0. 015 units smaller due to the pion excess. For instance, 
it reaches c 2 s = 0.127 at E lab = 40A GeV. Both mod- 
els indicate that at the energy around Ei ab = 4QA GeV 
the slope of the c 2 (y/s) distribution is changing, and the 
velocity of sound becomes less sensitive to rising bom- 
barding energy. 

Let us discuss the obtained values of the c 2 . For the 
ultrarelativistic gas of light particles the well-known the- 
oretical result is c s — 1/\M of the speed of light [28| . 
As shown in [29| , the presence of resonances in particle 
spectrum generates the decrease of the sonic speed. Em- 
ploying the empirical dependence [30| 

p(m) cx m a ' (2 < ol < 3) , (17) 

where p(m) dm denotes the number of resonances with 
masses from ratom + dm, one arrives to the equation of 
state in the form [29[ 

e=(a' + 4)P, (18) 

1. e., j < c 2 < g. This trend is reproduced in microscopic 
models. 

Since neither energy density nor pressure can be di- 
rectly measured in the central area of heavy-ion colli- 
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FIG. 10: (Color online) Time evolution of the microscopic 
pressure P and the energy density s in the central 125 mi 3 cell 
in (a) UrQMD and (b) QGSM calculations of central Au+Au 
collisions at energies from 11.6j4 GeV to 158^4 GeV. Dashed 
lines correspond to the non-equilibrium stage of the reaction, 
solid lines represent the equilibrium phase. 

sions, the experimental evaluation of the c s is more dif- 
ficult. One may rely on the hydrodynamic calculations, 
which claim that the magnitude of the so-called ellip- 
tic flow V2 depends on the speed of sound c s [31]. Us- 
ing the estimates, obtained for fixed impact parameter 
b = 8 fin under assumption of constant c s throughout 
the system expansion, PHENIX collaboration reported 
the value c s w 0.35 ± 0.05 [H| , i.e. c 2 s w 0.12 ± 0.3, for 
gold-gold collisions at top RHIC energy = 200^4 GeV. 
This value is close to our results and also implies rather 
soft effective EOS. 

Lattice calculations [33| predict an asymptotic value of 
c 2 s ~ 0.3 slightly below the Stefan-Boltzmann limit which 
indicates the appearance of a strongly coupled partonic 
medium. Recombination processes decrease the mean 
free path of the particles and lower thus the viscosity 
of the medium. By including such processes the sonic 
speed can be increased above the critical energy density 
thus coming closer to the lattice predictions [161 ]. 

The velocity of sound defines the change of entropy and 
energy densities with decreasing temperature, provided 
the local equilibrium is maintained during the expansion. 
The analytic expressions, which can be derived, e.g., for 
gas of non-strange mesons with zero chemical potential, 



FIG. 11: (Color online) The same as Fig. [10] but for macro- 
scopic pressure P extracted from the SM fit to microscopic 
data. 

read (see Appendix |Bj 

< 19 > 

The ratios e/eq and s/sq as functions of T/Tq obtained 
from model calculations at Ei a b = 20A GeV and 40^4 
GeV are plotted in Fig.[T3]together with results for fi = 
given by Eqs. (fl9|) - ((20| . Although the hadron gas in the 
cell represents a cocktail of species with different chemical 
potentials, that can be either zero, positive or negative in 
case of antiparticles, the curves calculated by the UrQMD 
and QGSM are not far from the ideal ones. Moreover, 
there is just a very weak difference between the UrQMD 
and QGSM curves for both energies. If one formally fits 
these distributions to Eqs. (fTT)l) - (f2"0"]) using the velocity of 
sound as fitting parameter, one gets a = 0.2 exactly. It 
would be nice to check whether the deceleration of en- 
ergy(entropy) density falloff with dropping temperature 
could be charged solely to the presence of hadrons with 
non-zero chemical potential. One way to do this is to 
perform a similar analysis of the cell conditions at RHIC 
(or higher) energies. - Here strange hadrons, baryons and 
their resonances are still present [ll| , but both chemical 
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FIG. 12: (Color online) The ratio P/e = c 2 3 , where P is 
defined microscopically, in the central cell of volume V = 
125 fm 3 as a function of center of mass energy %/s in UrQMD 
(solid symbols) and in QGSM (open symbols) calculations. 
Lines are drawn to guide the eye. 



potentials, /ie and /is, are quite small. Therefore, one 
may expect that the microscopic results would be closer 
to those presented by Eqs. ([T( | -([20] ) . 

Note also, that pressure in the cell changes with energy 
density quite smoothly, and no peculiarities which can be 
attributed to first-order phase transition are seen in the 
early stage of the reaction. - Here we simply extend the 
formalism of extraction of the thermodynamic parame- 
ters to the non-equilibrium phase, where one cannot trust 
the obtained values anymore. This was done merely in 
order to find any traces of the transition related to the 
onset of equilibrium and to changes of the effective EOS 
in the models. However, the analysis is performed for 
the fixed cubic cell of relatively large volume V = 125 
fm 3 , where the matter is distributed non-homogeneously 
at early times. To get rid of the evident ambiguities, the 
scheme is properly modified. 



D. Early stage of the evolution 

The central cell is further subdivided into the smaller 
ones, embedded one into another. The size of the initial 
test- volume is just Vmit — 0.125 fm 3 , and the energy den- 
sity e of the cells becomes the main parameter now. If the 
e of the inner cell is not the same (within the 5% limit of 
accuracy) as the energy density of the outer one, the SM 
analysis of the thermodynamic conditions is performed 
for the inner cell. If the energy density is uniformly dis- 
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FIG. 13: (Color online) (a) Ratio of energy densities e/eo vs. 
T/To in the central V = 125 fm 3 cell. Dashed line and solid 
line represent UrQMD calculations of central Au+Au colli- 
sions at 20A GeV and 40^4 GeV, respectively, whereas dotted 
line and dash-dotted line show the QGSM results for these 
reactions. Asterisks depict the analytic calculations given by 
Eqs. ([T9 ]l -p0 )l with fx = and a = c 2 s = 0.14. (b) The same 
as (a) but for the ratio of entropy densities s/so- 



tributed within the outer cell, the latter becomes a new 
test-volume, and so on. In the latter case it appears (see 
Fig. [9]) that the onset of the isentropic expansion regime 
in the central area occurs significantly earlier than the 
formation of equilibrated matter. Moreover, at the col- 
lision energies below 80 A GeV entropy per baryon ratio 
seems to be quite stable almost from the beginning of the 
reaction. 

Evolutions of the temperature and baryon chemical 
potential both in the central cell of the fixed volume 
V = 125 fm 3 and in the expanding energy area are de- 
picted in Fig. [T3J One sees 

that the transition to equilibrium proceeds quite 
smoothly if the analysis is performed for the fixed cell 
(Fig. 1 141 upper plot). In contrast, in the area with uni- 
formly distributed energy the transition to the equili- 
brated phase is characterized by a kink distinctly seen 
in each of the phase diagrams in both microscopic mod- 
els. Although this effect takes place along the lines of 
the constant entropy per baryon, it should not be auto- 
matically linked to the highly anticipated quark-hadron 
phase transition. The reason is simple, - extraction of 
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FIG. 14: (Color online) The evolution of the temperature T 
and baryon chemical potential [ib in the central cell of central 
Au+Au collisions at energies from 11.6^4 GeV to 158A GeV. 
Both parameters are extracted from the fit to the SM. Sym- 
bols and dashed lines show the evolution of these quantities 
in a cell of instantly increasing volume {Vimt = 0.125 fm 3 ), 
while dash-dotted (upper plot) and full (both plots) lines are 
related to calculations with the fixed volume V = 125 fm 3 . 



the thermodynamic parameters, such as T, P, /xb and /is 
(but not the entropy density, which is determined micro- 
scopically) , by means of the equilibrium statistical model 
is doubtful for the non-equilibrium phase. On the other 
hand, the formation of the kink may not be accidental. 
It is correlated with the significant reduction of the num- 
ber of processes going via the formation and fragmenta- 
tion of strings, and, therefore, with the inelastic (chem- 
ical) freeze-out of particles. - In both models the mat- 
ter, produced in a central area in central heavy-ion colli- 
sions at energies between AGS and SPS, is dominate d by 
(pseudo)elastic collisions after t w 6-^8 fm/c [33, l3ol l36j| . 
- In the fixed-cell analysis all parameters within the cell 
are averaged and the transition is smeared out. The ob- 
served phenomenon can easily mimic the signature of the 
QCD phase transition in the T-/iB plane, found in lattice 
QCD calculations (3?j] also along the lines of the constant 
entropy per baryon. 

Figure [T4l demonstrates also that thermodynamic char- 
acteristics of the fixed-size cell and the instantly growing 
energy-homogeneous area coincide completely during the 
equilibrium stage. In accord with earlier observation [Icl |. 



neither the mechanical reduction of the test-volume in 
longitudinal direction nor the criterion of uniformly dis- 
tributed energy density alone can help us in searching 
for quick equilibration in the central zone of relativistic 
heavy-ion collisions. Criteria of local thermal and chemi- 
cal equilibrium described in Sect. II VI are fulfilled after the 
chemical freeze-out in the test-volume, when the produc- 
tion of new particles in the system is ceased. 



VI. CONCLUSIONS 

In summary, two different microscopic string models 
were used to study the formation and evolution of the 
locally equilibrated matter in the central zone of heavy- 
ion collisions at energies from 11.6^4 GeV to 160^4 GeV. 
Calculations were performed both for the cubic central 
cell of fixed volume V = 125 fm 3 and for the instantly 
expanding area of homogeneous energy density. Tradi- 
tional approach based on the fulfillment of the conditions 
of kinetic, thermal and chemical equilibrium has been ap- 
plied to decide whether or not the equilibrium is reached. 
Both models favor the formation of the equilibrated mat- 
ter for a period of about 10 fm/c. During this period the 
expansion of matter in the central cell proceeds isentrop- 
ically with constant entropy per baryon. The equation of 
state can be approximated by a simple linear dependence 
P = a(y/s)e, where the square of the speed of sound 
c 2 s = a(y/s) varies from 0.13 (AGS) to 0.15 (SPS) in the 
UrQMD calculations, and from 0.11 (AGS) to 0.14 (SPS) 
in the QGSM ones. In both models the rise of a(y/s) with 
energy is slowed down after £) a (, = 40^4 GeV and satu- 
rates at RHIC energies. This change is assigned to the 
transition from baryon-dominated to meson-dominated 
matter. 

Study of the expanding area of the isotropically dis- 
tributed energy reveals that the relaxation to equilibrium 
in this dynamic region proceeds at the same rate as in 
the case of the fixed-size cell. However, the entropy per 
baryon ratio becomes constant before the state of equi- 
librium is attained. Here both microscopic models unam- 
biguously show the presence of a kink in the T-fiB phase 
diagrams. The higher the collision energy, the earlier the 
kink formation. Its origin is linked to the freeze-out of 
inelastic reactions in the considered area. 
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APPENDIX A: REDUCTION OF 
TEMPERATURE 

Let us consider non-relativistic ideal hadron gas which 
contains non-equilibrium number of pions, while the 
other hadron species correspond to their equilibrium val- 
ues. In thermal equilibrium the total energy of the gas 
is a sum of the masses of all particles (potential term) 
and the energies of their thermal motion (kinetic term). 
Compared to the case of fully equilibrized hadron gas, the 
temperature of the system with overpopulated amount of 
pions should reduce so that the total energies of both sys- 
tems remain the same. 

One can write 



(Al) 



= J2m^NP + 3 -T^J2 N ^ > ( A2 ) 



= ]Tmfivf + ^)^ N p, (A3) 



J, 7T "-"IT ' 

= N (2 J 



(A4) 
(A5) 



where the superscripts (1) and (2) are related to partially 
non-balanced (w.r.t. pions) and fully equilibrized system, 
respectively. Parameter a > 1 measures the excess of 
pions in system (1). From Eqs. (|Al|) - (|A5j) we have 



(T ( 2 ) _ T (i) J- flf) - ( a - l)T^N^ 

(A6) 



Introducing the reduced variables 

p = E^ 2) /M 2) , iP>i) 

i 

7 = , (0< 7 < 1) 

2 m w 

o 



3T( 2 ) 



we get finally 



/3(l- 7 ) = (a-l)(7 + S) 



(A7) 

(A8) 
(A9) 



(A10) 



Now knowing the pion abundance in particle spectrum 
(3~ l at chemical equilibrium and pion excess a one can 
estimate the drop of temperature 7 in the system due to 
re-distribution of kinetic energy among the extra degrees 
of freedom. 

APPENDIX B: EVOLUTION OF e AND s WITH T 

AT ,u = 

Gibbs free energy G is linked to energy E and entropy 
S of the system with pressure P, volume V and temper- 
ature T via the equality 



G = E + PV -TS 



(Bl) 



On the other hand, G = [xN, where \x is the chemical po- 
tential and N is the number of particles. If the chemical 
potential is absent, Eq. (|Blj) is reduced to the following 
expression for the energy and entropy densities, e = E/V 
and s — S/V , respectively: 



e + P = Ts 



(B2) 



Utilizing the condition li = 0. one can derive from basic 
thermodynamic equalities [28| 



de = T ds , 
dP = sdT . 



(B3) 
(B4) 



Inserting the equation of state dP = ade into these 
equations, we get after straightforward calculations 



a de dT 



and, finally, 



1+a e T ' 
ds dT 



- - (- 

so \T 

- - i- 



(B5) 
(B6) 



(B7) 
(B8) 



The obtained results are general for particles with fi = 
and do not depend on the expansion, e.g. longitudinal 
or spherical, scenario. 
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